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' We investigate the magnetic properties of the Cu-0 planes in stoichiomet- 

(N ' 

ric Sr„_iCu„+i02n (n=3,5,7,...) which consist of CuO double chains period- 
ically intergrown within the Cu02 planes. The double chains break up the 
two-dimensional antiferromagnetic planes into Heisenberg spin ladders with 



a 



' Ur = \{n — 1) rungs and ni = ^(n + 1) legs and described by the usual 



o 



antiferromagnetic coupling J inside each ladder and a weak and frustrated 
inter ladder coupling J'. The resulting lattice is a new two-dimensional trellis 

><" 

' lattice. We first examine the spin excitation spectra of isolated quasi one 

\ 

dimensional Heisenberg ladders which exhibit a gapless spectra when is 
even and ni is odd ( corresponding to n=5,9,...) and a gapped spectra when 
rir is odd and n/ is even (corresponding to n=3,7,...). We use the bond oper- 
ator representation of quantum S = ^ spins in a mean field treatment with 
self-energy corrections and obtain a spin gap of ~ ^ J for the simplest single 
rung ladder (n=3), in agreement with numerical estimates. We also present 
results of the dynamical structure factor S(q,a'). The spin gap decreases con- 
siderably on increasing the width of the ladders. For a double ladder with 
four legs and three rungs (n=7) we obtain a spin gap of only O.IJ. However, 
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a frustrated coupling, such as that of a trellis lattice, introduced between 
the double ladders leads to an enhancement of the gap. Thus stoichiomet- 
ric Sr„_iCu„+i02n compounds with n=3,7,ll,... will be frustrated quantum 
antiferromagnets with a quantum disordered or spin-liquid ground state. 
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I. INTRODUCTION 



The so-called "infinite layer" or "all layer" compound SrCu02 represents a new family 
of the Cu-0 superconductors. This compound crystallizing in the tetragonal structure is 
characterized by an infinite stacking of Cu02 planes with intervening Sr layers without 
oxygen representing the simplest possible charge reservoir. The infinite layer SrCu02 ^ are 
synthesized under extreme conditions involving high temperatures and high oxygen pressures 
and superconductivity in the Cu-0 layers is induced by modification of the intervening Sr 
layers. The Cu-0 layers can be either hole or electron doped depending on the dopants 
introduced in the Sr layer. The high pressure forms of SrCu02 are however unstable above 
certain temperatures and a homologous series of oxides formulated as Sr„_iCu„_|_i02n (with 
n=3,5,7,9, ) begins to be mixed in with the parent (n=oo) phase. 

The homologous series of Cu-rich high pressure phases Sr„_iCu„+i02„ were recently 
studied and were shown to consist of parallel lines of CuO-double chains periodically 
intergrown within the Cu02 sheets as illustrated in Fig. 1(a). The occurrence of these 
double chains was interpreted as resulting from a periodic shear operation with a shear 
vector of |(110) in the parent (n=cx)) phase. This can also be visualised as the appearance 
of domain walls within the planes. The structure of one such sheet for a general n is shown 
in Fig. 1(a). The Cu-atoms are shown as big black dots while the oxygen atoms are located 
at all the points of intersections of the straight lines. The Sr atoms (not shown in the figure) 
are located at the centers of the squares which are empty and in planes displaced by ± c/2 
with respect to the copper-oxide sheet shown with c being the lattice constant perpendicular 
to the Cu-0 planes. The dotted circles are the square co-ordinated Cu-atoms (like in the 
parent compound) and there are respectively 0,1,2,... such Cu-atoms in between two double 
chains for n=3,5,7... Thus the introduction of the CuO double chains periodically in the 
matrix of square planar co-ordinated Cu02 produces a superlattice geometry with unit cell 
parameters (naxa). 

The double chains in Sr„_iCu„4.i02r!, compounds affect dramatically the magnetic prop- 
erties of the copper-oxide planes. In the stoichiometric compounds the Cu sites are singly 
occupied (in the hole notation) and the 0-sites are empty. This can be modelled by a 
S=l/2 Heisenberg model on the new lattice, the trellis lattice, formed from the Cu-sites. 
The exchange interaction between Cu-atoms which are both not located along the double 
chains is given by that of the bulk (n=oo) value of J. Along the double chains however, two 
Cu-ions (such as A and B of Fig. 1(a)) are connected via an 0-site by 90° bonds which 
gives rise to a ferromagnetic exchange 0. The superexchange path between these Cu-ions 
is through two orthogonal 0-orbitals and this introduces a Hund's rule contribution leading 
to a ferromagnetic exchange J' (< 0) given in perturbation theory by 



A2 'Et + 2A Es + 2A 
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where tpd is the nearest neighbor hopping in the plane between the Cu3d(x^-y^) to the 
02p(x,y) orbitals, A = ep — with ep and being respectively the on-site energies of the 
02p(x,y) and the Cu3d(x^-y^) levels, and E^- and are respectively the triplet and singlet 
levels. Using the values of Ej- = 1.8eV and E^ = 7.3eV obtained from the level splittings 
1^ assuming an onsite Coulomb repulsion of Up =4eV and taking the standard values of 
tpd=1.3eV and A = 3.3eV [|[ we obtain JVJ=0.1-0.2. 

Thus the planes in Sr„_iCu„+i02n are broken up into Heisenberg ladders with a weak 
and frustrated interladder coupling as illustrated in Fig. 1(b). There are 0,1,2... vertical 
chains (separated by the width a) cutting the dashed region in Fig. 1(b) corresponding to 
n=3,5,7... . The lattice that results is a trellis lattice consisting of individual ladders with 
Ur = |(n— 1) rungs and n; = |(?t,+1) legs and coupled to each other through zig-zag couplings 
J'. Since J' <^J we can to a first approximation neglect J' and divide the Cu-0 planes into 

independent sets of ladders with 71^=1,2,3... rungs and n;=2,3,4 legs corresponding to 

n=3,5,7... . The compounds with odd n,, and even n^ (corresponding to n=3,7,ll,...) and 
those with even n^ and odd n^ (corresponding to n=5,9,13,...) exhibit different spin excitation 
spectra, the former having a spin gap with a shortrange exponentially decaying magnetic 
correlation function while the latter are gapless with a longrange powerlaw decay of the 
correlation function (like in the case of a single chain corresponding to n=l). 

It was argued in ref. p that Sr„_iCun+i02n compounds with n=3,7,ll... would be 
frustrated quantum antiferromagnets and spin liquids. This was motivated by the numerical 
studies 0, on isolated Heisenberg single- rung antiferromagnetic ladders with couplings J 
which exhibited a gap of ~ | J. Here we will follow an alternate treatment in terms of bond 
operators formulated to study two-dimensional dimerised Heisenberg systems. In section 
II we will emphasize some of the important features of these bond operaters and use them in 
a mean-field treatment to study the simplest single rung Heisenberg ladder with two legs. In 
section III we report the results of the calculations on the spin excitation spectrum, spin gap, 
ground state energy and the dynamical structure factor S{k,u>). We extend the calculations 
to double ladders and also to periodic ladders in section IV. We show that the spin gap 
decreases considerably on increasing the width (rungs) of the ladders. In section V we study 
the effect of frustration on two spin-ladders and show that any non-zero frustrated coupling 
leads to an enhancement of the gap. This important result points towards the stabilization 
of a spin-liquid ground state in a trellis lattice. Such lattices will be realizations of short 
range RVB (Resonating Valence Bond) ground states in a S=| system [0. 



II. SINGLE SPIN-LADDER 

We investigate here the properties of a single spin- | ladder, shown in Fig. 2(a) which is 
described by the standard Heisenberg antiferromagnetic model 
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H = J^Si^ ■ Sr^ + \J Sm^-Sm,+i, (2.1) 

i i,m=l,r 

where we take J to be the strength of the interaction along the rungs(i) of the ladder and 
A J the interaction along the legs of the ladder. S^^ and S^^ are respectively spin-i operators 
at the left(l) and right (r) -hand sites on each rung i of the ladder. In the limit of A=0 
(strong- coupling limit) the Hamiltonian of Eq.(2.1) reduces to a sum over contributions 
from independent two spin rungs. Thus it would be natural to tackle this problem from the 
limit of singlet dimers placed on the rungs and then switch on the interaction between them. 
It will be seen later that this starting point leads to reasonable results even for the case of 
A=l which is of interest here. We follow the bond-operator representation of quantum S=| 
- spins introduced to study specifically the properties of dimerized phases. We emphasize 
here the essential features of the bond operators. 

We consider two S = | spins Si and placed on each rung. The Hilbert space is spanned 
by four states which can be combined to form the singlet \s > and the three triplet \tx >, 
\ty >, and \tz > -states defined as being created out of the vacuum |0 > by the singlet and 
triplet creation operators 



\s > 


= s^\0 >= 


^(ITi>-UT>), 




\tx > 


= 4lo>= 


-^(ITT>-Ui>), 




\ty > 


= 4lo >= 


^(ITT>+I ii>), 


(2.2) 


\tz > 


= 4|o>= 


^(in>+UT>). 





A representation of the spins S/ and Sr in terms of these singlet and triplet operators is 
given by, 

Sia = ^{sha + tls - ieafB^tlt^), (2.3) 

Sra = 2^ ^^^a ^^0/374^7)' (2-4) 

where a, /3, and 7 represent respectively the components along the x,y, and z-axes and e 
is the Levi-Civita symbol representing the totally antisymmetric tensor. Henceforth, it is 
assumed that all repeated indices over a, /3, and 7 are summed over. 
A constraint of the form, 

sU + tit^ = l, (2.5) 

is introduced for each dimer in order to restrict the physical states to either singlets or 
triplets. Taking the singlet and triplet operators at each site to satisfy the bosonic commu- 
tation relations, 
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= 1, [t^,tl]=5^p, [5,4] =0, (2.6) 
we can reproduce the S=|, SU(2) algebra of the spins S/ and : 

Si-Sr = -ls^s + h%, Sf = S'^ = l. (2.7) 

Substituting the operator representation of spins defined in Eqs.(2.3) and (2.4) into the 
original Hamiltonian of Eq.(2.1) we obtain the following form: 



H = Ho + X{Hi + H2), (2.8) 



where, 



^0 = E ■Ji-l4s^ + \tlu^) - E t^Msi + tlu^ - 1), (2.9) 

i i 

Hi = +^J2 4ji+iasl+iSi + tljl^^^SiSi+i + h.c), (2.10) 

H2 = E 2 '^■'^ ~ '^a/3)(4a4+la^j/3^i+l/3 - 4a4+l/3^j+la^i/3 + h.c). (2.11) 

i 

The part of the Hamiltonian containing triple "t" operators vanishes identically in the 
present case due to reflection symmetry [jlT|. A site-dependent chemical potential fii is 
introduced to impose the constraint of Eq.(2.6). The Hamiltonian of Eq.(2.8) can now be 
solved by a mean-field decoupling of the quartic terms. This yields an effective Hamiltonian 
HmR with only quadratic operators. We take < Si >= s, which means that the "s" bosons 
are condensed. We replace the local constraint fii by a global one /i in accordance with the 
translational invariance of the problem along the ladder-axis (y-axis here) . We will consider 
here only the terms Ho and Hi in Eq.(2.8) and will show later in the Appendix that inclu- 
sion of H2 changes the results only slightly. We perform a Fourier transformation of the 
operators 4a = J2k Aa^^^^'^ where N is the number of dimers or rungs in the ladder and k 
is the wave-vector whose only non-zero component is along the ladder axis. Thus, retaining 
only terms Hq and Hi in Eq.(2.8) we obtain the mean-field limit, 

Hmifi, S) = N{-jjf - lif + /X) + Y\^ktljka + Afc(tL^lfca + hat-ka)l (2.12) 
^ k 

where, 

J 



Afc = — — yu + XJs^ COS fc, 
A 



/\k='-Js'^cosk. (2.13) 
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We have taken the lattice constant to be unity. The ground state wave-function can be writ- 
ten in the form |(^ >= Cexp[J2i Sis\ — Y.k ^fc^Ia^^fco] |0 >• The above mean-field Hamiltonian 
Eq.(2.12) can be diagonalised by a Bogoliubov transformation into new Boson operators 
given by, 

Ika = cosh Oktka + sinh 9kt^_f,^, (2. 14) 

where the coefficients cosh^^ and sinh^^ obtained in terms of A^, and ujk are given by, 

cosh^^,^ 1(^ + 1), 

sm\i^ek = \sgn{^k){—-l)- (2.15) 

We finally obtain, 

Hm{^^, s) = N{--^Js^ - fxs^ + fx)- ^(^ - f,) + J2M^L^ka + I), (2.16) 

k 

where, 

LJk=[A\-{2Aky]K (2.17) 
The parameters fj, and s are determined by solving the saddle point equations; 

<^>=0,<^>=0. (2.18) 

OfX OS 

We obtain the following self-consistent eqations from Eq.(2.18) evaluated at T=0. 



where K(?7) and E{r]) are respectively the complete elliptic integrals of the first and second 
kind of modulus 77. The dimensionless parameter d is defined as 

d-jrEiry (2-20) 



III. RESULTS 
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A. Spin-triplet spectrum 



For each A, the self-consistent solutions of s and /i are obtained from Eqs.(2.19) and are 
used to determine the excitation spectrum of the system of a single ladder obtained from 
Eq.(2.17) as, 

= J(l _ /f)[l + (icosfc]5. (3.1) 

These quasiparticle excitations arising from the spin-triplet states of the spin-ladder form a 
band whose bandwidth is a function of A. The band-minimum is at k=7r and in Fig. 2(b) we 
present plots of uJk/J (relative to the band minima) as a function of k for A=0. 1,0.5 and 1.0. 
We notice that the dispersion around the band minimum gets more linear with increasing A 
which is reminiscent of the case of a linear chain (A = oo). 

Recently extensive numerical calculations have been performed []T2| on Heisenberg spin 
ladders using Lanczos techniques. The spin-triplet dispersion relations were obtained on a 
ladder with 2x12 sites and for various values of A. The results of such a calculation are 
shown as filled circles in Fig. 2(b) and the agreement with the present work is very good. 
For very small values of A the spin-triplet excitation spectrum has a bandwidth of 2AJ which 
is in excellent agreement with those obtained from the Lanczos data and from the strong 
coupling expansions of Heisenberg ladders [0. The dispersion relation of Eq.(3.1) can be 
parametrized by a spin-wave velocity which is given by Cg = J{j — j)(f)^- The spin-wave 
velocity reduces to as A — >0 and is I.IJ for A = 1.0. 



B. Spin Gap 

If the excitation spectrum Uk is real and positive everywhere in the Brillouin-zone then 
the system is in a magnetically disordered (spin-liquid) phase and has a spin gap given by, 

A = J(l-^)[l-rf]i (3.2) 

This is indeed found to be the case for a single Heisenberg antiferromagnetic ladder. In 
Fig. 2(c) we present a plot (continuous line) of the spin gap (in units of J) as a function 
of A obtained numerically by solving the non-linear Eqs.(2.19) and then substituting the 
self-consistent solutions s and into Eq.(3.2). Alternatively, Eq(2.19) and Eq.(2.20) can be 
combined to give the following single equation for d 




which is easier to solve numerically. The solution of d for a given A can be used in Eq.(2.19) 
to determine fi and the spin gap is then obtained from Eq.(3.2) 
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We can analytically study the asymptotic behavior of the spin gap. For small values of 
A (and hence d) the elliptic integrals of Eq.(2.19) can be expanded in a power series as, 

V/; 2^ 4' 64' ' 
K{v) = l{l + -/ + -■■■). (3.4) 



and we obtain, 



d = 2\[l-l\^ + 0{\% (3.5) 

o 

j) = l + ^A2 + 0(A^), (3.6) 

and the spin gap is given by 

A = J(l-A-^ + ^ + 0(A^)). (3.7) 

This is consistent with Fig. 2(c) which shows a linear drop of the spin gap A for small 
values of A. As A gets larger the spin gap as shown in Fig. 2(c) deviates considerably from 
the linear behavior and does not exhibit any critical value for A where the gap vanishes. 
Thus starting from the limit of A = with spin singlets on each rung of the ladder, we see 
that any finite coupling A along the legs of the ladder delocalizes the singlets on the rungs 
thereby reducing the magnitude of the gap but not closing it completely. 

The spin gap seems to approach zero as A ^ 00 which should be the right gapless limit 
corresponding to decoupled spin-| chains. However, on closer inspection of Eq.(3.3) we 



notice that at A = 00, d is given by the solution of the equation -7=K(y^^) — ^ = 0. 
This gives a value of d which is very close to 1. In addition it can be shown using Eq.(2.19) 
that as A ^1, j diverges as Aln(l — d). Hence the spin gap given by Eq.(3.2) also diverges 
as A — >• cxD. In the present case the spin gaps start to increase for values of A > 3 and the 
mean-field treatment ceases to be valid. 

The Heisenberg ladder has also been extensively studied numerically P], [0] using the 
Lanczos technique and the spin gaps have been determined on 2xN ladders with N=4,6,8,10, 
and 12. An extrapolation of these results to the bulk limit is shown (as filled circles) in the 
inset of Fig. 2(c) as a plot of versus A~^ for values of A < 1. The full curve in the same 
inset indicates results obtained from the present calculation. The agreement is found to be 
good only for small values of A. 

The deviation of the present spin gaps from those obtained numerically can be traced to 
the A^ terms in the expansion of the spin gap as obtained in Eq.(3.7). In Ref. |]12| it was 
shown in a strong coupling expansion that the spin gap varied as A = J(l — A + |A^) for 
small values of A. The A^ term arises in the strong coupling expansion from short range 
effects as the two nearest neighbor singlets surrounding the rung which is excited to a triplet 
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state are ineffective in contributing to the energy of spin-triplet excitations to O(A^). Thus 
the positive co-efficient obtained for the A^-term shifts the spin gaps above the 1-A line as 
A is increased in agreement with the numerical results (see filled circles in the inset of Fig. 
(2c)). However, in the present treatment a positive coefficient is obtained in the expansion 
for (i — ^) (Eq.3.6) but the expansion of (l-d)^ produces an overall negative co-efficient for 
the A^-terms which shifts the spin gaps below the 1-A line as A is increased. This implies 
that for a finite A the singlet and triplet levels are not pushed apart enough due to the 
wrong treatment of the short-range effect in the present mean-field method. This arises 
from treating the local constraint to be valid only on the average. We have therefore added 
an additional self energy term of O(A^) in the triplet levels of Eq.(2.9) which corrects for 
the neglect of short range effect in the present treatment. We use the self enrgy term (3)? 
with an optimum value of /3 = 0.7 which gives reasonable values of the spin gaps for A < 1 
as shown by the dashed line in the inset of Fig. 2(c). It should be pointed out that this self 
energy correction does not modify the dispersion curves of the spin-triplets shown in Fig. 
2(a), it merely shifts the positions of the minima at A; = tt. 



C. Ground State Energy 

The ground state energy obtained from Eq.(2.16) is given by (neglecting H2 term) 



§ = - + « - i(i - « + i(i - m + ^M^^j). (3-8) 

where d is defined in Eq.(2.20) and /i = j- 

In the limit of A ^ the correct ground state energy ^ = — | is recovered. As A 
is switched on the ground state energy decreases and we present in Table I the values of 
2^ obtained from Eq.(3.8) for certain values of A. As discusssed before for the spin gap 



the ground state energies have also been obtained numerically [|T2| on finite length ladders 
and the extrapolated values to the bulk limit is also presented in Table I. We observe that 
the energies obtained from the present mean-field treatment compare well with numerical 
estimates for A < 0.5. 



D. Structure Factor 

Structure factors are important physical quantities which can be measured experimen- 
tally. The dynamic spin-structure factor is defined as 



S^'^'ig^u) = J dte'-' < S'^'J{t)S':lM >, (3-9) 
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where S^'^^ represents the Fourier transform of the cu-th component of the spin-operator 
combination along the rungs given by 

S'at) = j:^''^-^'[SUt) ± Sr^t)], (3.10) 

i 

where S/. and S^^ are the two spin-operators at each rung i of the ladder as shown in Fig. 
2(a). Using the representations of the spins S; and in terms of the singlet and triplet 
operators as defined in Eqs.(2.3)and (2.4) we obtain 

Substituting these in Eq.(3.10) and transforming the "t" operators to Bogoliubov operators 
defined in Eq.(2.14) we obtain the following expressions for the dynamic spin structure 
factors; 

S"^(g,a;) = \f{cosh2eq - sinh2eg){ng + e{uj))5{ujg - \uj\), (3.12) 

S"{q,uj) = ^ ^{[cosh2{9k+q - Ok) + l]nfe(l + nk+q)S{uJk+q - ujk - u;) 
^ k 

1 

+ -[cosh2{9k+q - 9k) - l](nfc + Q{uj)){nk+q + Q{uj))5{uJk+q + i^k- kl)}, (3.13) 

where,©(a;) is a step function and n^ is the Bose occupation factor and cosh^^ and sinh^fe 
are defined in Eq.(2.14) 

The first term of Eq.(3.13) represents the simultaneous emission and absorption of ex- 
citations and vanishes identically at T=0 while the second term corresponds to creation 
or annihilation of two excitations. In Fig. 3 we present the results of the structure factors 
S^'^{q, uj) as a function of oj for q=7r and for A = 1.0 after applying the self-energy correction 
described at the end of section 111(b). The dominant contribution to the structure factor 
comes from S^{q,Lu). The static structure factors Sgl^{q) = J dujS^'^{q,uj) show peaks at 
the commensurate wavcvcctor of q= n. The spin-spin correlations in the present case decay 
exponentially at large distances with a correlation length given by, 

f = l = t2(T^l*- 

where Cg is the spin- wave velocity and A the spin gap and d is as defined in Eq.(2.20). We 
notice that ^ — > as A — > and for A = 1 we obtain a value of ^ = 2.66. 



IV. DOUBLE AND PERIODIC ARRAY OF LADDERS 

We now consider the effects on the spin gap of increasing the number of rungs in the 
ladder. We first consider the case of a simply connected double rung ladder shown in Fig. 
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4(a). The two spin-ladders (denoted by left(L) and right(R)) are connected with a strength 
of A' J and in each individual ladder we assume as before an interaction strength of J along 
the rungs and AJ along the legs of the ladders. The Hamiltonian of the system can be 
written as, 

H^T.[iHo, + XH,^) + lx'H[^], (4.1) 
where, i denotes the ladder index, Left(L) or Right(R), and 

~ 2" 't\iat\+Ua^H^i+U + h.C.)^ (4.2) 

i 

Here, we have neglected terms of the form H2 of Eq.(2.11) as it is shown in the appendix 
not to change the results significantly in the parameter range < A < 1, which is of interest 
here. It should be pointed out here that the Hamiltonian containing triple "t" operators 
will not vanish here in the H' part of the Hamiltonian connecting the two ladders as it did in 
the previous case along the ladder axis. However this will not give any contribution in the 
magnetically disordered phase when we take averages over the operators in the mean-field 
decoupling. 

Replacing < >= s and ^jLi^ — ^ va. Eq.(4.2) and taking Fourier transforms of the 
t-operators we obtain 

H^ili, s) = 2N{-- Jf - ,if + ^) + Y.[^ktljk,a + ^k{t\J_k,^ 

+tkeat-kea) + ^{'^^kea'^ke/a + tkea'^-ke/a + tkeat-ke/a)]- (4.3) 

where and are as defined in Eq. (2.13) and = — ^s^J. We perform a Bogoliubov 
transformation into two new Boson operators defined in terms of the t-operators of the left 
and right hand ladders as 

7i,2,„ = ^ [(cosh Oktk^^ + sinh Oktlf^^J ± (cosh Okiki^a + sinh ^jttLfc^^)] . (4.4) 

These are simply symmetric (bonding) and anti-symmetric (anti-bonding) combinations of 
the individual transformations in the left and right ladders. Wc can now diagonalise the 
Hamiltonian of Eq.(4.3) using this transformation and we obtain the following 

Hmifi, s) = 2N{-Ijs' - fif + ^^)-Nd-|,)+ ^ M^fl^^,,^ + h, (4.5) 

^ ^ k,b=l,2 ^ 
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where 0Ji^2k is defined as, 

001,2, = [Afc' - (2A,)2 ± 2Tk{Ak - 2Afc)]i (4.6) 

Thus the spin-triplet excitation spectrum of a double ladder with three rungs and four legs 
consists of two branches for each A and A' representing the bonding and the anti-bonding 
states. The splitting of these two branches is governed by the transfer matrix proportional 
to Ffe which is in turn proportional to A'. As A' — > the two branches collapse into a single 
branch of the one rung ladder described in section Ill(a). We have plotted in Fig. 4(b) 
the spin-triplet excitation spectrum for the double ladder with the vahic A = A' = 1 and 
after applying the self-energy correction in each ladder as described at the end of section 
Ill(b). For comparison we show in the same figure a plot (shown as the dashed curve) of 
the excitation spectrum of a single rung ladder (A' = 0). 

Following the same procedure as described in section II we can write down the mean-field 
equations evaluated at T=0 as 

2 _ 3, _ 1 /■ dky l + ia+ l + jg- . 



where d is defined in Eq.(2.20), and a± is given by 



a±^cosk±—. (4.8) 
4A 

The expression for the spin gap is 

A = J(^-^)[l-4l + ^)]i (4.9) 

We have plotted in Fig. 4(c) the spin-gaps of the double ladder as a function of A' for 
different values of A. We obtain a value of O.IJ for the spin gap for A = A' = 1 and we notice 
that for each A the spin gap reduces drastically with A'. This can be explained by the fact 
the singlets along the rungs of the two ladders can not only delocahze along the ladder axes 
via the A-coupling but also across the ladders through the transfer matrix proportional to A'. 
This produces for any non-zero coupling of two single rung ladders (A' 7^ 0) a splitting in the 
excitation spectrum of the individual ladders into two branches, one above (the antibonding 
branch) and the other below (the bonding branch) as shown in Fig. 4(b). Since the spin 
gap is a measure of the minimum energy in the excitation spectrum (at k = it) and since 
one of the branches of the double ladder is always lower than that of a single ladder for any 
non-zero A' we can conclude that the spin gap of a double ladder is lower than that of single 
ladder and the spin gap should progressively decrease on increasing the width of the ladder. 
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In the same context it is worthwhile to consider a periodic array of ladders,with A=l in 
each ladder, and simply connected to one another by A' and then to ask the question for 
what value of A' should the spin gap disappear since it is known that the ground state of 
a two-dimensional antiferromagnct has longrange order. We can study the system of lad- 
ders arranged periodically within the formalism described previously and the self-consistent 
equations are given by 

2 3, 1 f f d^k 1 + la 
''--2^ = -2 J J 



(2x)2 [l+da]i 

'2 



- + 2- =+A//- r^r r- 4-10 

where k is now a two-dimensional wave-vector with components ky along the ladder-axis 
and kx across the ladders. The parameter d is define in Eq.(2.20) and a is given by 

a = cos A;^^ — — cos /Cj;. (4-11) 

Combining the equations in Eq.(4.10) we can write down the following single mean field 
equation for d which is easier to solve numerically. 

d f f d^k 1 , . ^ 

- = 3.0-//-— r r- 4.12) 

The excitation spectrum of the periodic system of ladders has a minimum at k = (tt, tt) 
and the the spin gap given by A = J(| — — d{l + ^)]^ is plotted in Fig. 5 as a function 
of A' with A set to 1.0. We have again applied the self-energy corrections on each individual 
ladders as described at the end of section Ill(b). We notice that the spin gap vanishes for a 
value of A' 0.25. As described before for the double ladder the decrease in the spin gap 
with A' is explained by the delocalization of the singlets across the ladders. In this case the 
transfer matrix connecting the ladders is twice as large as in the case of a double ladder 
thereby decreasing the spin gap with A' even faster than that of a double ladder. 



V. FRUSTRATED DOUBLE LADDERS AND TRELLIS LATTICE 

We now consider the double ladder described in the previous section but connected to 
each other through zig-zag frustrated couplings J' = X' J as shown in Fig. 1(b), where J is 
the standard rung coupling. We retain the Hamiltonian in the form of Eq.(4.2) but with 
H']^^ containing an additional term given by 

~ X] i^ia^i+Wot^k^i+'i-il + AiaA-^Ula^k^i+ltl + h.c), (5.1) 
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The mean field Hamiltonian and the excitation spectrum are those described by Eqs.(4.5) 
and (4.6) with Tk replaced by a new form = — J(l+cos k) We obtain the self-consistent 
equations of Eq.(4.7) but with a-t given by a± = cos k ± ^{1 + cosk) 

The excitation spectrum of the double ladder with a zig-zig frustrated coupling also 
consists of two branches like in the case of a simply connected double ladder. However, 
in the present case at the minimum position of the spectrum {k = n) the two branches 
become degenerate. This can be explained by the fact that the singlets on two successive 
rungs of each ladder are completely out of phase and hence are not able to delocalise across 
the ladders through the zig-zag frustrated couplings A'. Thus the spin gap defined by the 
minimum value of the spectrum at /c = vr is now given by, 

A = J(^-^)[l-rf]^ (5.2) 

which does not contain the A' term explicitly. Thus the spin gaps should not change much 
form the values obtained for a single ladder and this is indeed what is observed as shown 
by the dashed lines in Fig. 4(c). We have again applied the self-energy correction to each 
individual ladder. We notice from Fig. 4(c) that for each A the spin gap slightly increases 
with A'. This simply reflects the slight changes in the self-consistent solutions of d and /j, 
with A'. 

We can extend the results discussed above to the trellis lattice shown in Fig. 1(b). The 
excitation spectrum would now be a function of k with components ky along the ladder axis 
and k^ across the ladders. At the point ky — ir the spectrum would be dispersionless with k^ 
as a result of destructive interference of successive rung-singlets along the ladder axis. Thus 
a frustrated coupling between the ladders in the trellis lattice will not affect the spin gaps 
of each individual ladder and hence it will only help in retaining the spin-liquid nature. 

VI. CONCLUSIONS 

We have shown that the CuO double chains intergrown periodically in the Cu02 planes 
of the new series of infinite-layer compounds Sr„_iCu„+i02n (n=3,5,7,...) creates a new two- 
dimensional spin-lattice, the treUis lattice. Such a lattice can be described by Heisenberg spin 
ladders with = |(n — 1) rungs and ni — + legs with the usual antiferromagnetic 
coupling J inside each ladder and a weak and frustrated inter ladder coupling J' <^J. On 
neglecting J' the Cu-0 planes can be thought of as built up of independent quasi one- 
dimensional ladders with odd n^ and even n; (corresponding to n=3,7,ll,...) and those with 
even n^ and odd n/ (corresponding to n=5,9,13,...) which exhibit different spin excitation 
spectra. The former are gapped with shortrange exponentially decaying magnetic correlation 
functions while the latter are gapless with a longrange powerlaw decay of the correlation 
functions. 
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We have used the bond operator representation of quantum spins in a mean field treat- 
ment with self-energy correction to first study the excitations of the simplest single rung 
ladder with two legs corresponding to n=3. We have obtained the spin-triplet dispersion 
with a minimum at = tt and a spin-gap of (~ These are in good agreement with 

numerical estimates 0, |jl2[. It should be pointed out that conventional spin- wave calcula- 



tions applied to a ladder ||T3[, [T^ predict gapless excitations for all values of the couplings 
A J along the legs of the ladder, even for A=0 which is clearly unphysical. 

We have been able to extend the mean field treatment to double ladders and to a periodic 
array of ladders. We find that increasing the rungs of the ladder drastically reduces the spin 
gap. For a double ladder with four legs and three rungs corresponding to n=7 we obtain a 
spin-gap of only O.IJ. In a periodic array of ladders with intraladder coupling J we find that 
the spin gap vanishes for an interladder coupling of ~ 0.25J. 

We have also studied the effect of a frustrated coupling, such as that of a trellis lattice, 
introduced between two ladders. We find a slight enhancement in the spin gap. Extending 
the results to a trellis lattice we can show that the spin-liquid nature will be preserved. Thus 
stoichiometric Sr„_iCu„_|_i02n compounds with n=3,7,ll ... will be frustrated quantum 
antiferromagnets with a quantum disordered or spin-liquid ground state. The resulting 
trellis lattice will be a realization of the short range RVB ground state for a spin-i system 



The implications of these results to the other high Tc compound can be conjectured 
as follows. Certain underdoped high Tc samples have been experimentally shown to have 
spin gaps and a theoretical description of this in terms of frustrated next nearest neigbour 
coupling in a two-dimensional antiferromagnetic lattice is obtained only with large and un- 
physical values of the coupling. However if the two dimensional Cu-0 planes were thought of 
having some microstructure (introduced upon doping) then the spin gaps could be explained 
quite naturally. 
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APPENDIX: 

The results presented in section III were obtained by neglecting the term H2 of Eq. (2.11) 
containing four triplet operators. Here we will study the changes on including this term in 



16 



a similar mean field treatment as presented before. By taking quadratic decouplings of the 
operators in H2 and performing Fourier transformations we obtain, 

k 

where A and N are as defined previously and P and Q are two new mean-fields defined as 

P — < tloti+la >, 

Q = < tiati+la > ■ (^2) 

On including Eq.(Al) into Eq.(2.11) we obtain a mean field Hamiltonian similar to 
Eq.(2.12) 

+ Hi^ktijka + MiJ-ka + tkat-ka)], (A3) 
k 

but with Afe and replaced by 

J 2A 
Ak — — — /J, + XJs^ cos k — 3"-^*^ '^^^ ^' 

Ajfc = +-Js^cos A; + ^QJcos A;. (A4) 
2 3 

We can diagonalisc Eq.(A3) by a Bogoliubov transformation as described earlier and the 
spectrum is given by Eq.(2.17) but with and A^ replaced by their new forms. The 
parameters 11, s, P and Q are obtained by solving the saddle-point equations which reduce 
at T=0 to the following equations; 

dk Afe 



^ 2^ J 2Ti2ojk 
/3^„/u f dk Ak-2Ak , 



i^k 

dkA^ 
2% 2cuk 

Q = f f^— cos/c. 

J ZTT OJk 

By numerically solving for the mean-field parameters the spin gap is obtained and is plotted 
in Fig. 2(b) (dashed lines) as a function of A. We notice that the inclusion of H2 does not 
change the results significantly even for A f=:i 1 
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FIGURES 

FIG. 1. (a) Schematic diagram of a single copper-oxide sheet in Sr„_iCu„+i02n showing the 
parallel lines of CuO double chains. The Cu-atoms are shown as big black dots while the oxy- 
gen atoms are located at all the points of intersections of the straight lines. The dashed regions 
correspond to the usual square co-ordinated Cu02 regions, (b) The two-dimensional trellis lattice 
formed from the exchange couplings in a single copper-oxide plane of Sr„_iCuri_|_i02n at stoichiom- 
etry. The usual antiferromagnetic coupling J describes each of the ladders having n^ = ^(n-1) rungs 
and n; = ^(n-|-l) legs and the ladders are coupled to one another through weak and frustrated 
zig-zag couplings J'. 

FIG. 2. (a) A single S=^ Heisenberg antiferromagnetic ladder with couplings J along the rungs 
(i) and AJ along the left(l) and right(r) legs of the ladder, (b) The dispersions of the spin-triplet 
excited states of the ladder of (a) relative to the band minimum at k=7r for several values of A. 
The continuous curves are obtained from the present mean-field method and the filled circles are 
the Lanczos results of Ref.l2 obtained on a 2 x 12 ladder, (c) The spin gap A (in units of J) as 
a function of A obtained from the present mean-field treatment (continuous curve). The dashed 
curve is a similar plot after including higher order terms as described in the Appendix. In the inset 
a plot of A/A J versus A~^ is shown. Continuous curve is from the present treatment (without 
higher order terms) and the filled circles are the numerical results of Ref.8. The dashed line in the 
inset is obtained after including a self-energy correction [described at the end of section 111(b)] to 
the present treatment. 

FIG. 3. The dynamical structure factors S'^''^{q,(jj) as a function of frequency lo for q=7r. 
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FIG. 4. (a) A double Heisenberg ladder consisting of two single S=^ antiferromagnetic ladders 
of Fig. 2(a) connected to each other through the couplings A' J. (b) The dispersions of the 
spin-triplet excited states (bonding and anti-bonding states of Eq.(4.6)) of the double ladder of (a) 
for A = A' = 1.0. The dashed curve is a similar plot for a single ladder (A' = 0) with A = 1. (c) 
The spin gap A (in units of J) as a function of A' for several values of A. The dashed curves are 
similar plots obtained for the case of a double ladder connected to each other through the zig-zag 
frustrated couplings J' = \' J of Fig. 1(b). 

FIG. 5. The spin gap (in units of the rung coupling J) of a periodic arrangement of ladders as 
a function of the inter ladder coupling A' and with an intraladder couling of A = 1. 
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TABLES 



TABLE I. Results of the ground state energy 21^7 of a Heisenberg ladder 



A 




Eg 
2N.J 








Present work 




Ref. 


12] 


0.5 


-0.394 




-0.43 


1.0 


-0.475 




-0.578 
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